Band Structure and Quantum Conductance of Nanostructures from 
Maximally-Localized Wannier Functions: 
The Case of Functionalized Carbon Nanotubes 
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We have combined large-scale, F-point electronic-structure calculations with the maximally- 
localized Wannier functions approach to calculate efficiently the band structure and the quantum 
conductance of complex systems containing thousands of atoms while maintaining full first-principles 
accuracy. We have applied this approach to study covalent functionalizations in metallic single- 
walled carbon nanotubes. We find that the band structure around the Fermi energy is much less 
dependent on the chemical nature of the ligands than on the sp^ functionalization pattern disrupt- 
ing the conjugation network. Common aryl functionalizations are more stable when paired with 
saturating hydrogens; even when paired, they still act as strong scattering centers that degrade the 
ballistic conductance of the nanotubes already at low degrees of coverage. 
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Carbon nanotubes (CNTs) are some of the most 
promising materials in micro- and nano-electronics, 
whose applications could range from transistors to inter- 
connects to memories 1]. In the present effort to engi- 
neer and tune the electronic properties of CNTs, organic 
functionalizations of the sidewalls provide one of the most 
encouraging avenues 0, Q , ideally combining the excep- 
tional electronic properties of CNTs with the diversity of 
organic chemistry. In this Letter, we study the effects of 
covalent functionalizations on single-walled CNTs using 
a novel approach to calculate the electronic structure and 
quantum conductance, based on our idea of constructing 
well-localized and transferable orbitals as ideal building 
blocks for the electronic structure of complex nanostruc- 
tures 0, 0, S 0- These maximally- localized Wannier 
functions (MLWFs) provide a most compact and most 
accurate local representation of the electronic structure 
of a solid or a molecule, and are used to construct the 
full Hamiltonian and to derive the band structure and the 
Landauer conductance of complex systems containing at 
least thousands of atoms, preserving first-principles ac- 
curacy and with linear-scaling computational costs. 
Pristine, single-walled CNTs are ideal one-dimensional 
quantum wires with quasi-ballistic electron transport for 
tens of nanometers (at low bias) even at room tempera- 
ture 1^ . Structural or electronic defects and inelastic 
excitations are then some of the most relevant factors 
which limit conductance. Covalent functionalizations 
will alter conduction manifolds and doping levels, thus 
affecting significantly the electronic transport properties 
of CNTs [lOj. First-principles approaches are especially 
suited to characterize such effects, in the controlled en- 
vironment of quantum-mechanical simulations. Here, we 
use density-functional theory in the GGA-PBE approxi- 
mation using state-of-the-art solid-state approaches 



. Explicit first-principles calculations of systems con- 
taining thousands of atoms are rarely possible 13], and 
the delocalized nature of the Bloch orbitals is not well- 
suited to describe the transport properties of aperiodic 
systems; these obstacles are overcome by the approach 
described in the following. 

Disentanglement and Localization: Our method starts 
from the construction of an optimal MLWFs basis set. In 
a large supercell, where the sampling of the Brillouin zone 
(BZ) is limited to the F-point, MLWFs can be obtained 
by a unitary transformation of the occupied ground-state 
eigenfunctions |w„) = Y.ni=i ^rnn iV'm), where 

U is chosen to minimize the sum of the spread of all 
around their centers 0, Q|. This is not a viable ap- 
proach for metallic systems, or even for insulators where 
the occupied subspace can significantly change charac- 
ter as we move away from F, even if for the supercell 
at hand F-sampling might be adequate to calculate the 
energy and charge density of the system with the nec- 
essary accuracy. Our goal is instead to obtain, from a 
combination of the occupied and some of the unoccu- 
pied eigenfunctions at F, a compact and localized set 
of orbitals that describes a group of bands of interest 
(i.e., manifold) in the whole BZ of the supercell - sim- 
ilar in spirit to the k • p method but with a dif- 
ferent approach and improved outcome. This optimal 
set is obtained with a two-step procedure. First, an 
optimally-smooth subspace is extracted from the Hilbert 
space (i.e., "disentangled"), using a procedure recently 
introduced and specifically adapted here to the case 
of F-sampling. This iV-dimensional subspace consists of 
a linear combination of all No eigenstates at F that fall 
inside a given energy window (with No > N) , and is cho- 
sen so that its dispersion across the BZ is as small as 
possible: this is akin to say that its projection operator 
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FIG. 1: MLWFs obtained from the disentanglement and local- 
ization procedure. We show in (a) the 3 inequivalent MLWFs 
obtained for a pristine metallic (5,5) CNT, clearly correspond- 
ing to a chemical picture of pz and sp^ orbitals, and in (b) 
the same MLWFs for the case when an aminophenyl group 
(1% coverage) has been covalently attached. 

has a minimal change of character as the BZ is swept, 
or that the basis set obtained is as complete as possi- 
ble. This subspace is identified by its projection operator 

P^^^ = J2n=l\M{(t>n\, where |(/)„) = Em°=l ^™"IV'm)- 

The No X N matrix A (such a.sA'^A = I) is the key 
quantity to be determined, and is chosen so that the pro- 
jection operator at F has maximal overlap with its first 
neighbors, i.e., it maximizes Y,i lFiTr[P(i')p('^')], where 
the neighbors G; and the weight factors Wi are chosen 
according to Refs. 0,^3- Once this optimally-smooth 
manifold is extracted, a standard localization procedure 
can be applied on the {ji/in}} to obtain MLWFs {\wn)}. 
Band Structure and Quantum Conductance: The pro- 
cedure outlined above corresponds to an inverse map- 
ping of first-principles, extended Bloch orbitals into tight- 
binding orbitals |w,ir) = |w„(r — R)) that are then peri- 
odically repeated to generate, via Bloch sums, wavefunc- 
tions at any arbitrary k: iV'nk) = 7^ Z]r e*'^'^|w„R) 
(R and k refer to the direct and reciprocal spaces of the 
supercell). Localization allows us to neglect in the re- 
sulting Hamiltonian matrix all terms {wmR\H\wnii') for 
which |r„iR - r„R/|>L, with r^R the center of \wmii), 
and L a cutoff distance determined by the spread of the 
MLWFs. Since these MLWFs span both the occupied and 
unoccupied subspaces, they will be well-localized and de- 
cay exponentially even in metals. The Hamiltonian 
matrix becomes diagonally-dominant, and is straightfor- 
wardly given by 

(V'mkl-fflV'nk) = {WmQ\H\Wno) + ^ 6*"^'^ (WmO | ^ | W„r) 

R 

(1) 

where the sum runs only over the very few supercells 
(two, for 1-dimensional nanotubes) interacting with the 
one at the origin. The advanced and retarded Green 
functions G^'^ of the conductor and the couplings with 
the left or right lead T^^ji can be derived from the ma- 
trix elements of the Hamiltonian, providing the equilib- 
rium Landauer conductance C/(£') = 2e_Xr(rLGpF_RGp) 
(see Refs. H E H [13 for detailed explanations of 
this approach and its implementation in periodic bound- 
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FIG. 2: Band structure, quantum conductance, and density 
of states of a (5,5) CNT as calculated in the present approach 
based on MLWFs (solid lines) , and compared with the results 
of a full diagonalization (circles). The band structure has 
been unfolded over the BZ corresponding to the primitive 20- 
atom unit cell; the five vertical dashed lines indicate the 5 
k-points corresponding to F in a 100-atom supercell. 

ary conditions). As a validation, we first calculate the 
band structure and the quantum conductance of a pris- 
tine (5,5) CNT in a 100-atom supercell. The manifold 
of interest is that spanned by the bonding combination 
of sp^ orbitals (the graphitic backbone) and by the Pz 
orbitals (tt and tt* bands), so the disentanglement and 
localization procedure uses a target dimension of 2.5 or- 
bitals per atom. We show in Fig. the resulting ML- 
WFs, confirming this intuitive chemical picture. In ad- 
dition, the right panel shows the case of an aminophenyl 
functionalization, where the "half-filled" Pz orbital is re- 
placed by an sp^ bonding orbital. The two other ML- 
WFs barely show any change, clearly hinting at the fact 
that maximal-localization corresponds also to maximal- 
transferability of the MLWFs building blocks. Once the 
(short-ranged) matrix elements of the Hamiltonian have 
been determined, the band structure at any point in the 
BZ can be obtained with negligible computational costs 
by diagonalization of Eq. ^ . Results are shown in Fig. 
121 and compared with the band structure obtained from 
a diagonalization in a complete plane- wave basis. The 
agreement between the two calculations is excellent - es- 
sentially a proof of the proper disentanglement of the two 
manifolds (represented by solid lines and circles) and lo- 
calization of the MLWFs. The extrema in the band struc- 
ture and the steps and peaks in the quantum conductance 
and the density of states are also all in perfect agreement. 
Even more significantly, the massless bands, the quan- 
tum conductance, and the density of states around the 
Fermi energy clearly display the metallic character of the 
nanotube, notwithstanding the fact that the MLWFs have 
been obtained from the eigenstates at F in the presence of 
a pseudo-gap of eV in the folded BZ. The reason for 
this agreement is that, moving away from F, the change 
in character of the occupied manifold is described exactly, 
for all practical purposes, by the unoccupied part at F of 



our extracted manifold. 

These tools now allow us to study with full first-principles 
accuracy very complex structures, where the MLWFs ob- 
tained from a calculation at T or even from an isolated 
fragment (e.g., a hydrogenated nanotube fragment) will 
reproduce accurately the fine details of a much larger ex- 
tended system. We focus here on metallic CNTs whose 
sidewalls are functionalized for tens of nanometers with a 
disperse array of covalent ligands. Functionalization with 
aryl groups has been recently accomplished [2flj , and our 
goal is to characterize the effects of electronegative or 
electropositive ligands such as nitro- or amino-phenyls 
on the band structure and transport properties of CNTs. 
It is important to stress that the addition of a covalent 
ligand can affect nearby sidewall carbons, and expedite 
further reactions that reduce the local strain energy or 
minimize the disruption of the conjugated manifold, as 
observed in fuUerenes lUl- While exhaustive studies are 
necessary to fully characterize the thermodynamic stabil- 
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ity of the functional group decorations [22 
subtle spin polarization they might induce 
that, once a functional group (FGi) is covalently bound 
to the sidewall, the attachment of a second one (FG2) to a 
nearby carbon becomes more facile. These reactions are 
strongly site-dependent; we present in Fig. |31 the first 
and second attachment energies, AEi and AE2, corre- 
sponding to CNT + (FGi)H ^ CNT(FGi) -I- I/2H2 and 
CNT(FGi) + (FG2)H ^ CNT(FGi)(FG2) + I/2H2 re- 
spectively 2Ql . In order to understand how the electronic 
structure of the tube is modified by the presence of the 
ligands, we first studied a periodic array of functionalized 
pairs in the configuration (a) of Fig. 13 at a density of 1 
pair per 100 carbons. The resulting band structures are 
shown in Fig. ^ for the cases where the functional group 
pairs are composed either by two hydrogen atoms, or a 
nitro- or aminophenyl and a hydrogen. For comparison, 
we also report the band structure of the pristine CNT 
(folded in the BZ of the supercell), and the one obtained 
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FIG. 3: AE2 in eV as a function of the position of the func- 
tional group FG2 when the group FGi is attached in the po- 
sition marked by X (A_Ei is indicated in parenthesis) on a 
(5,5) CNT. FGi and FG2 are, respectively, (a) H and H, (b) 
nitrophenyl and H, (c) aminophenyl and H, (d) phenyl and 
phenyl (in our calculation, the nitro and amino residues are 
in the para position). The arrow lies parallel to the nanotube 
axis, and we denoted with a, 13, 7, and 5 the arrangements 
for four of the most stable pairings. 




FIG. 4; Band structure of a (5,5) CNT decorated with a peri- 
odic array of functionalization pairs in position a (1 pair every 
100 carbons): (a) H pair (dots are from a Quantum-Espresso 
calculation p^. for comparison), (b) nitrophenyl and H, (c) 
aminophenyl and H, (d) model calculation, where the 2 pz 
MLWFs on the functionalized sidewall carbons have been re- 
moved from the calculation, (e) pristine nanotube. 



from the pristine CNT but removing the two Pz MLWFs 
sitting on the two sidewall carbon atoms at (a) . Interest- 
ingly, these band structures are largely insensitive to the 
chemical nature of the attached groups, notwithstanding 
the charge transfer taking place when e.g., aminophenyl 
is substituted for nitrophenyl (0.1 |e| per 100 atoms), or 
the chemical differences between a phenyl moiety and a 
hydrogen. The foremost effects on the bands around the 
Fermi energy are described by a simple model in which 
the functionalized sidewall carbons change hybridization 
from sp^ to sp'^ . Functionalization transforms the Pz 
MLWF into a bonding sp^ MLWF whose on-site energy, 
roughly lowered by 7 eV, removes it altogether from the 
TT conduction manifold; this effect is clearly captured by 
the model (Fig. Eld)). A perfectly periodic array of func- 
tional pairs might be quite unrealistic; from the trans- 
port point of view, the fundamental quantity of interest 
is the reduction in the transmission probability as more 
and more groups are attached to a functionalized seg- 
ment of an otherwise pristine nanotube (see the top of 
Fig. 0). Our approach is particularly suited to study 
such complex cases. We show first in the left panels of 
Fig. 13 the Landauer conductance as a function of energy 
for an infinite (5,5) CNT with either one pair attached 
in positions (a), (/3), (7), and {5) or a single covalent 
functional group. The covalently-bonded groups act as 
strong scattering centers, reducing the conductance at 
the Fermi energy by 19%, 8%, 21%, 17% and 42% re- 
spectively. Again, the dominant factors are geometric 
and similar results are obtained, say, from hydrogen at- 
tachments (dashed line) or the removal of Pz MLWFs on 
the relevant sidewall carbons (solid line) . Once a random 
array of functional groups covers the nanotube for tens of 
nanometers - comparable to the lengths of pristine tubes 
for which we can observe ballistic transport |^ - the con- 
ductance of the tube drops dramatically, as shown in the 
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FIG. 5: Top: infinite metallic (5,5) CNT functionalized in its 
central region by an array of phenyl pairs. Left panel: quan- 
tum conductance of an infinite (5,5) CNT with one isolated 
pair of ligands in positions a (a), /3 (b), 7 (c), 5 (d), or with 
a single ligand attached (e) (dashed line: ligands are hydro- 
gen atoms; solid line: model calculation; dotted line: pristine 
CNT). Right panel: random distributions of pairs of ligands 
or of single ligands (arranged as in the left panel), for the 
case of 10 defects per 1000 carbons (solid line) or 30 per 3000 
(dashed line), averaged over 5 random configurations. 



right panels of Fig. |S1 for the case of a functionahzed 
conductor 12.3 or 37.0 nm long (corresponding to 1000 
or 3000 carbons) between semi-infinite leads. We note 
that substitutional dopants such as boron and nitrogen 
also introduce substantial (albeit weaker) disorder, and 
can change the transport regime from ballistic to local- 
ized li^. Among the patterns studied, (3 turns out to 
be the weakest scattering center, since it is a paired de- 
fect that preserves the mirror plane containing the tube 
axis and is almost transparent to either of the two Bloch 
states with well-defined parity over the mirror plane (a 
similar result has been discussed for symmetric Stone- 
Wales defects on a (10,10) CNT H^). 
In conclusion, we have presented a novel approach to 
calculate the band structure and transport properties 
of complex nanostructures from large-scale electronic- 
structure calculations at the F point via an essen- 
tially exact inverse mapping into explicit and localized 
tight-binding orbitals (MLWFs). These MLWFs act as 
electronic-structure building blocks, representing at the 
same time a minimal and most accurate basis set for the 
Hilbert space at hand and providing a natural avenue 
to map first-principles calculations into model Hamilto- 
nians. An application of this approach to covalently- 
functionalized metallic CNTs provides a clear correla- 
tion between chemical bonding and transport properties, 
showing that the band structure around the Fermi en- 
ergy is strongly influenced by the sp^ functionalization 
pattern disrupting the conjugation network. Common 



aryl functionalizations were found to be more stable when 
clustered with saturating hydrogens; even when paired, 
they still act as strong scattering centers that degrade 
the ballistic conductance of carbon nanotubes already at 
low degrees of functionalization. 
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